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Abstract 

A self-interacting polymer can undergo an orientational ordering 
transition, depending on the magnitude of the nematic interaction. 
The effect of embedding such a polymer into a flexible surface on this 
transition is studied on the mean-field level. Renormalized values of 
the elastic constants of the "dressed" surface are derived as functions 
of the orientational order parameter of the polymer chain. In the 
disordered state the surface tension and curvature modulus remain 
scalars but depend on the surface coverage of the embedded polymer. 
In the nematic state there is a symmetry breaking transition leading 
to anisotropic elastic constants. At a sufficiently large nematic order 
parameter the effective surface tension in the direction perpendicular 
to the nematic axis can become negative, leading to tubulization of 
the embedding surface. 
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1 Introduction 



Usually the adsorption of polymers onto surfaces is discussed in terms of 
more or less diffuse adsorption layers, composed of trains, loops and dan- 
gling ends. There is enough corroborating evidence that in many cases this 
simplified picture correctly captures the salient features of the polymer ad- 
sorption phenomena 

However, some recent results on adsorption of highly charged polyelec- 
trolyes onto oppositely charged molecular surfaces give evidence to much 
stronger bonding of the adsorbing polymer to the surface [0, [|. It seems that 
the polymer adsorbs flat onto the surface without any observable loops, dan- 
gling ends and even a very restricted number of crossings among the different 
segments of the chain, as the strong electrostatic attraction can apparently 
override any entropy driven tendency for less tightly bound configurations. 

In this respect one can envision the adsorption process more as an em- 
bedding of polymer chains into (or onto) the adsorbing (supporting) surface. 
The theoretical significance of this observation would be that an analysis of 
the adsorption phenomenon could be carried out without taking into con- 
sideration the extremely improbable configurations with local desorption of 
the polymers (loops and dangling ends). In certain sense this simplifies the 
analysis because the chains can be treated as embedded into the supporting 
matrix. However, the embedding also signifies that there is a strong coupling 
between the polymer and the surface configurations and thus the pertaining 
degrees of freedom can not be decoupled and treated separately. 

In this contribution we shall address the interplay of the polymer and sur- 
face degrees of freedom in connection with the nematic ordering transition 
of a polymer chain embedded (i.e. tightly adsorbed) onto a flexible, fluid 
surface. It will be established that there can be no effective decoupling of 
the tight adsorption problem (embedding) into a separate polymer and sur- 
face parts. The influence of the surface degrees of freedom on the polymer 
ordering transition, as well as of the polymer degrees of freedom on the effec- 
tive elastic constants of the "dressed" (meaning surface + polymer) surface 
will be derived and discussed. The renormalized form of the surface elastic 
constants will be shown to lead to a tubulization transition of the embedding 
surface. 
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2 Summary 



I first present a sketchy overview of the main results and conclusions reached 
in the paper. Starting from an orientationally (nematically) interacting poly- 
mer chain embedded into a soft (membrane) surface, conditions for a second 
order symmetry breaking orientational ordering transition are derived. The 
presence of surface shape fluctuations only marginally influences the critical 
point of the transition. The nematic orientational order parameter, defined in 
a standard way || as S = \ (3<cos 2 6> — 1), where is the angle between 
the local direction and the nematic axis, is obtained as a solution of 



S 



\ V 2(/3ep)/X 



(1) 



where e is the elastic modulus of the polymer chain, p is the polymer surface 
density (surface coverage) and u' a is proportional to the excluded volume of 
the nematic interaction u a , renormalized due to the presence of the embed- 
ding surface shape fluctuations, 

(3u' a — ► f3u a - \-^<{ci - c k f>Q. (2) 

where q is the embedding surface curvature in the direction i, where i = x,y. 
The average in the above equation is with respect to the bare surface, i.e. 
with no embedded polymers. 



Tied to this transition are also changes in the overall ("dressed") elastic 
constants of the membrane. A closed expression is derived for the renormal- 
ized elastic constants, surface tension (7) and the elastic modulus (K c ), of 
the "dressed" membrane 

7 — ► 7 + 2-^<7 Q + ec a 
K c — > K c + ^a a a/3, (3) 

where the index a describes the components of the indexed quantities in the 
reference frame of their eigen-coordinate system. Thus a can be either n 
signifying a component parallel to the nematic axis, or _L signifying a com- 
ponent perpendicular to the nematic axis. a a is the mean-field orientational 
tensor of the embedded polymer, and c a its mean-field squared curvature 
tensor. A is a scalar depending on the orientational order parameter. 



Following from the above expression for the renormalized values of the 
elastic constants, we show that the renormalized surface tension can change 
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sign in the direction perpendicular to the nematic axis for a sufficiently large 
nematic order parameter. Basing our analysis on the recent general investi- 
gation of the tubulization transition of fluid membranes by Radzihovsky and 
Toner || we show, in the framework of a Landau theory, that the preferred 
conformation of the surface is a tubule with the long axis perpendicular to 
the nematic direction of the embedded polymers. The tubulization transition 
of the membrane does not coincide with the isotropic - nematic transition of 
the embedded polymers. 



3 Analysis 



3.1 Outline 



The formalities of the derivation are a bit convoluted so I will give a short 
guided tour before actually developing them. The point of departure is the 
usual flexible chain Hamiltonian together with the "embedding ansatz" , con- 
straining the chain to lie on the supporting surface. The expression for the 
partition function following from these premises is developed in Section 3.2. 
A closed form of the partition function can be obtained by introducing collec- 
tive variables of the polymer chain: density, orientational tensor and tensor 
of the squared curvature. This is developed in Section 3.3 together with the 
final form of the partition function in collective coordinates in Section 3.4. 
This partition function has a simple closed form solution if the collective 
variables are treated in a mean-field manner, ignoring local fluctuations in 
their values along the surface. This is done in Section 3.5 leading to the final 
derivation of effective elastic constants of the "dressed" membrane in Section 
3.6. Once these are derived the tubulization transition can be described by 
means of the formalism set forth by Radzihovsky and Toner |5[ as is done in 
Section 3.7. 



3.2 Theory and Formalism 



The position of the n — th bead along the polymer contour is described by a 
position vector r{n) = (x(n),y(n), z(n)), which is a continuous, differentiable 
function of the arclength n, while the total length of the polymer is N. I 
limit myself to a single polymer chain case, but a generalization to many 
chains is quite straightforward. Let the embedding surface be represented 
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in the Monge parametrization as z = ((x,y) = C(p)> where p is the two- 
dimensional radius vector. In this notation the partition function of a flexible, 
but otherwise noninteracting, polymer chain embedded into the surface ((p) 
can be written as 

E (C (p(n))) =/-../ Vr{n) J] S (z(n) - ((p(n))) J[ 5 (f {nf - l) xe^ f ^\ 

J J n n 

(4) 

where 

W>(?(n)) = / n dn=^el o (f(n)) dn. (5) 



/3e is the energy associated with polymer bending and can be related to the 
persistence length 0. By integrating out the z(n) variable and enforcing the 
continuity condition r(n) 2 = 1 only globally [|7|, i.e. 



'r{nf) = {p{nf + am) 2 ) = ^ (6) 
one is led to the following form of the corresponding partition function, 

H (C (p)) = / • / Vp(n) exp (-(3H (p(n))) , (7) 

where some irrelevant constant terms have been omitted. The effective 
Hamiltonian Tlo(p(n)) has been obtained || as 

rN rN 

(3H (p{n)) = l(3ej Q (p*{n) 2 + ({n} 2 ) dn + Xj (p» 2 + ({n) 2 ) dn 

= \fr [ [?{nf + Hn) + |^| Mn)Mn)X) dn + 

v/. d((p)d((p) 



where the constant A is the Lagrange multiplyer, assuring the global condition 
Eq. |5], to be determined later. The validity of the above Hamiltonian is 
restricted to the case of a global continuity constraint Eq. ^ (for details of 
this transformation see f| and |§). 

In the case of a selfmteracting polymer chain the Hamiltonian has two 
additional terms corresponding to the interactions between different segments 

rN .. rN rN 

f3H (r(n)) — ► ±/3e / r(n) 2 dn + \(3 / dndm V iso (r(n) - rim)) + 

Jo Jo Jo 

nN . . - 2 

dndm (r(n) x f(m)J V an - iso (r(n) - r(m)) , 

(9) 



(8) 
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where the interaction has been split into an isotropic part and an anisotropic 
part, of a general nematic form [^]. 



After integrating out the variable z{n) and taking into account the fact that 
the chain is embedded into the surface C(p) one obtains for the interaction 
part 

dndm V iso (p(n) - p(m)) + ±/3 / / dndm T{p{n), p{m)) V an _ iso (p(n) - p{m)) 

Jo Jo 

where I used the shorthand V (p(n) — pirn)) = V {pin) — p{m); ({p{n)) — £(p(m))) 
and defined 

F(p(n), p{m)) = (p{n) 2 + ( 2 (p(n))) (p{mf + C 2 (p(m))) - (?{n)p{m) + C(p(n))C(p(m))) 2 = 

•2-2 

= Pi (n)p k (m) - p^p^p^njp^m) + 

/ • 2 • • 2 • • • • • \ d( d( 

+ ^Pi(7n)p fe (n)p,(n) + p (n)p fe (m)p(m) - 2p i {n)p i {m)p k {n)p l {m)J — — 

In the definition of J-{p{n), p{m) I limited myself only to terms up to and 
including the second order in the derivatives of CO)> m order to be con- 
sistent with the harmonic approximation assumed later for the bare surface 
Hamiltonian. 



+ . . . 
(11) 



The final integration in the partition function has to be over the surface 
degrees of freedom C(p). Assuming that in the absence of the embedded 
polymer the surface energy is given by V{({p)), the partition function can 
be cast into the form 

E(N) = J...J VC{p) S(C(p)) x exp{-f3V(C{p))). (12) 

In what follows I start from the assumption that V{({p)) is a quadratic 
functional of the form 

v(C(p)) = \ij (v ± C(p))Vp + \k c J (viC(p))Vp = 

i / d((p)d((p) r a 2 C(p) d^jp) 2 

= 2 / ^ik—^Z -^— d P+ 2 / K iklm ^ ^ r,^ r,^ » PA 13 ) 

J dpi dp k 2 J dp i dp k dp l dp m 

where 7 is the bare surface tension, and K c is the bare curvature elastic 
modulus. In order to facilitate later computations the surface energy and 
the curvature modulus tensors are introduced, which for a bare isotropic 
membrane assume the straightforward form ^ ik = 7^ and K ik[m = K c 5ikSi m . 
By gathering all the terms in In E(£(p)) that contribute to the quadratic order 
in C(P) one w iU t> e i n position to evaluate the rescaled values of the elastic 
modulus and the surface tension, which is the aim of this computation. 
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Also conveniently for later developements we introduce Sq(Q) through 
Fourier decomposed £(p) 

v(C(Q))=EMQ)\((Q)\ 2 , (14) 
Q 

with 

S Q {Q) = \iQ 2 + \K c Q\ (15) 

So(Q) is of course connected with the structure function of the bare mem- 
brane. 



3.3 Collective Variables 

One proceeds by defining collective variables in terms of which the partition 
function Eq.[7] assumes the form that allows an explicit introduction of mean 
- fields. The three collective variables introduced to this effect are: density, 
orientational tensor and the tensor of the squared curvature 

rN 

P(p) = / dn5(p-p(n)) 
Jo 

rN 

Vik(p) = / dn pi(n)p k (n) 6(p-p(n)) 
Jo 

rN 

Cik(p) = / dnp i {n)p k {n)8{p-p{n)). (16) 
Jo 



It should be noted here that since all the above collective variables are defined 
on a 2D surface, they do not share the same properties with their 3D coun- 
terparts. For instance the orientational tensor does not satisfy the equality 
Tr (Jik(p) = p{p) |0] • Later on an analogous relation valid in 2D will be 
derived in detail. Furthermore the squared curvature tensor is not connected 
with the torsion of the polymer chain but can be simplified by means of the 
Frenet equation for a planar curve into 



rN 

[p) = I dn K 2 (n) fii{n)n k {n) 8{p - p{n)) = </? 2 (n)n i (n)n fc (n)> c/iain , 
Jo 

(17) 

where n(n) is the (planar) curvature and Hi(n) is the normal vector of the 
chain. 



Cik( 



With collective variables the partition function Eq.|7| can be cast into the 
form 

S(C(p)) = /• • •/ Vp(n)Vp(p)Va ik (p)V Cik {p) x 



6 ~ L dn x 

<5 (vik{p) - J o dn ^(n)p fe (n) 5{p - p{n))^j x 

sUik(p)-J o dn p^njp^n) 5(p - p(n))j x 

x exp(-(3H(p(n),p(p),a lk (p),c ik {p))). (18) 



The Hamiltonian H (p(n) , p(p) , <Ti k (p) , Qfe(p)) is obtained by substituting 
the collective variables into 7io(p). Assuming furthermore (for the sake 
of convenience) that both the isotropic as well as the anisotropic parts of 
the interaction potential have a form of a contact interaction Vi SO (p, p) = 
UiS 2 (p — p) and V an _ iso (p, p) = u a 5 2 (p — p) one remains with the following 
form of the interaction Hamiltonian 

rN .. rN 

(3H(p(n),p(p),a tk (p),c ik (p)) = \(3e dn p(nf + A j dn p(n) 2 + 
+ &Jd P dB d , c lk (p) + Xjd p a ik (p) + 



dpi dp k J dpi dp 

+ \fr fd 2 p p-\p) f5ff p@ a ik {p)a lm {p) + 
J dpidp k dp l dp m 

+ (3u a I d 2 p[au(p)a k i(p) - cr ik (p)au(p)] ^——r + 
J dp k dp t 

+ j d 2 p p 2 (p) + \j3u a J d 2 p [a 2 u {p) - Tr (a lk (p)) 2 ] . (19) 

Here only terms up to and including the second order in ((p) and its deriva- 
tives have been retained. A higher order term can be derived from the 
Pi(n) p k {n) pi(n) term in Eq.|8| after the p{n) integration and has been conse- 
quently omitted. 



3.4 The Partition Function 



The next step is to introduce functional Fourier transform representations for 
delta functions appearing in Eq.[H| For each delta function we thus obtain an 
auxiliary field the integration over which insures that the definitions in Eq.|TB| 
are satisfied. Calling the auxiliary fields 0(p) for p(p), i>ik(p) for a ik (p) and 
ti k (p) for Ci k (p) we obtain the following relation 

S(CO)) = J ...J Vp{p)Va ik {p)Vc ik {p) x Vcj){p)V^ ik {p)Vt ik {p) x 
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Z(4>(p),i/Jik(p),tik(p)) x exp(-^(C(p);0(p),^fc(p),^fc(p);p(p),^fc(p),Q Jk (p) )), 

(20) 

where 

s {<f>{p), ipik{p), tik{p)) = /•••/ v p( n ) ex P OH; 0(0), ^(p), Ukip)))- 

(21) 

The definitions employed in the above two formulas were 

rN .. rN . rN 

^W o (p(n);0(p),Vifc(p),*ifc(p)) = |0e / dn p{nf + X dn p{nf - X dn + 

Jo Jo Jo 



rN rN . rN 

+ i dn <j>{p{n)) + i dn V>ifc(0(n))0i(n)0 fe (n) +i / dn £j fc (0(n)) 0^71)0^) 



jO JO 



(22) 



while 



^(CO);^^),^^),^^)^^),^^),^^) ) = 

ifl />* (X\ ^(P)^C(P) , . />-, m dC(p)dC(p) 
- 2 (3ejd p cM ^ ^ + X Jd p a lk (p) ^ ^ + 

5 2 C(p) #C(p) 



<9C 0C 



+ 

m 



+ 



+ |/9e Jd 2 p p 1 (p) a ik (p)o lm (p) 

+ (3u a d 2 p[a ii (p)a k i(p)-a ik (p)a i i(p n 

J dp k dp l 

+ | d 2 p p 2 (p) + |/3w a y d 2 p [a 2 (p) - Tr (a ifc (p)) 2 

- i Jd 2 p p{p)4>{p) - i Jd 2 p (T ik (p)ip ik (p) - i Jd 2 p c ik {p)t ik (p). 

(23) 



The partition function of the membrane with an embedded flexible polymer 
chain can now be obtained in several straightforward steps. Let us first of 
all introduce the renormalized potential V(C(0)) which is of the same form 
as Eq.[13] but with renormalized values of 7^ and K ik i m 

lik(p) — ► J$ik + tc ik (p) + 2^a ik (p) + 2u a [cr u (p)a ik (p) - cr K (0)cr Zfe (0)] 

Kikim{p) — > K c 5 ik 5 lm + -(p)a ik (p)(Ji m (p). 

P 

(24) 

The renormalized values of 7^ and K ik i m thus depend on the set p(0), cr ik (p) 
and Cjfc(p). The corresponding Fourier representation would ensue as 



V(C(Q)) = \Y,(lik{Q)QiQk + K iklm (Q)Q i Q k QiQ m j |C(Q)| 2 . (25) 

Q 
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Now we construct a generating functional (functional Fourier transform) of 
exp —ftV(C(p),p(p), <Tik(p), Cjfc(p)) at some quenched value of the membrane 
shape CO)- By definition one has 

exp (-/3V(C(p), (j){p),^ ik {p),t ik (p)fj = Z^fy ((f)(p),^ ik (p),t ik {p)) =/•••/ T>p(p)Va ik (p)Vc ik (p) x 
x exp ^-/3\/(C(p),p(p),a ife (p),c ife (p)) -i jd 2 p p(p)(f)(p) -i Jd 2 p a ik (p)^ ik (p) -% Jd 2 p c ik (p)t ik [ 



where the l.h.s. should be read symbolically as a shorthand for the functional 
Fourier transform. 

Finally the partition function can be expressed as an average of the prod- 
uct of the generating functional with variables 4>(p),ipi k (p),ti k (p) and the 
partition function of a single polymer chain in the orienting field of external 
sources 4>{p) , ipi k {p) and t ik (p), averaged over all the conformations of the 
membrane, i.e. 

E(N) = J VC(p)ZW)) (27) 

with 

Z(t(P)) =/••/ V(P(p)Vij ik (p)Vt ik (p) exp(-/3y(C(p),0(p),^fe(p),^(p))) E(<j>(p),i/; ik (p),t ik (p)). 

(28) 

Obviously the general form of the partition function is untractable and ad- 
ditional approximations have to be considered to get a closed-form solution. 



3.5 Mean-field Approximation 



It was recently shown p) that a similar partition function can be evaluated ex- 
plicitely by introducing a mean - field description of the polymer collective co- 
ordinates p(p), a ik (p), Cjfc(p) as well as their auxiliary fields 4>(p),ipik(p),ti k (p) 
by supressing their p dependence and treating them as constants. Thus 
by introducing jo(p) -> p,cr ik (p) -> a ik ,c ik (p) -» c ik as well as i(f)(p) -» 
0, iipik(p) — > ipi k , iti k (p) — > i ik one can first of all explicitely evaluate the 
partition function of a single chain with external sources H $i k , Uk) . This 
calculation will not be repeated as it has been already detailed in [[|. The 
final result valid in the limit N ^> 1 can be obtained as 



—kT In H (0, ^ ifc , Uk) = yJ2 



(29) 
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where index a refers to the eigenvalues so that ^ a ,i a are the indexed eigen- 
values of ipik and t ik . The next step now is to evaluate -Z^p) (^,^ik:Uk) on 
the mean - field level. In complete analogy with the calculations in || I 
obtain 



d 2 p +^o- a jd 2 p 



dp c 



+ 



+ \(3e p 1 a a (Tp Jd 2 p 



d 2 ((p) d 2 ((p) 



+ 



dp* m 

+ f3u a (Tra Q ) a p J d 2 p ( J - i5u a a a afi J d 2 p 



dCjp) gC(g) 
dp* dps 



+ \f3iiiS p 2 + \(3u a S (Tra. 



S <J a ip Q S c a t c 



2 

- s 



- n2 



Tr (a a ) 2 



(30) 



What remains now is the final integration over the membrane modes C{p) 
or equivalently over their Fourier components ((Q)- It is performed in such 
a way that all the terms describing the bare membrane and thus depending 
only on S (Q) are discarded as they do not enter the mean - field evaluation 
of collective variables or their auxiliary fields. The final expression obtained 
in this way for S(iV) thus assumes the form 



- kT\nE(N) 



= IE 

- s 



(A + ^a) 



\ (f + *c 

- S °o$* - S J2 c»t a + \l3uiS p 2 + \f3u a S (Tra a ) 2 - Tr (a a ) 



+ 



Q 



Qi 



PS Q (Q) 



+ 2AE 



Q: 



QlQl 



H — / ^ a ag- 

PS (Q) PtZ PSoiQ) 



(31) 



In the above equation I have limited myself to the lowest orders in c a and 
a a , meaning to the first order in Q 2 terms and to the second order in a a for 
Q 4 terms, all the higher orders have been discarded. This essentially limits 
the validity of conclusion derived from Eq.[H] to the regime of low surface 
coverage of the membrane, i.e. to small p. 



Introducing now 



(5e 



PS(Q) = (3S (Q) +PeJ2 cMl + 2A ]T a a Q\ + ~ E WpQlQp (32) 

a a P a , /3 



I obtain the following set of mean-field equations (the leftmost quantity 
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stands for the minimizing variable) 

Cry • Vr 



fo v Ql 

2 q PS(Q) 



N/S 



i> 



\ + ip a 

(3u a [(Tra a ) - c 



] + Aa a 

Q 



- > -1 N/S Q 



n (3e 

Ql 

(3S(Q) P , 



1 + 



2L 

(3e 



71 - V~* QlQ% 

+ — E^E 



Q 



(3S{Q) 



N/S 



^(X + lp*) (f + t a ) 



p = N/S 

^E- 



(X+^Pa) (f +t a ) 



a Q 



Q 



(3S{Q) 



(33) 



The appearance of wave- vector q max is a consequence of the limit of the con- 
tinuous description of a polymer chain and should be inversely proportional 
to the Kuhn's length. 



At this point a perturbation expansion in terms of the surface coverage of 
the polymers p is introduced and only the lowest order is taken into account, 
essentially making the approximation S(Q) — > S (Q). Since S (Q) is by 
assumption isotropic one has 

p^P \ dp a dp a / 
P^P \ dp a dp p I Q 

(34) 

Index signifies that the averages are taken with respect to a bare membrane. 
Obviously F is the strength of the surface normal fluctuations and F\\ and 
F± are the strengths of the surface curvature fluctuations. 



\ ~> Qa 

Q 



PS{Q) 



E 

Q 

E 

Q 



Ql 



PS (Q) 

QlQl 

(3S {Q) 



At this point the calculation performed here substantially deviates from 
M. Since in that case the polymer was ideal, i.e. non-interacting, the orien- 
tational ordering was due to higher order expansion terms in S(Q), omitted 
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here, as they would only rescale the already existing transition point but not 
its qualitative characteristics. 



The mean-field equations can now be solved by means of the ansatz 

^ = 27ITf) (1±s) ' (35) 

where S is the orientational order parameter S = \ (3<cos 2 0> — 1) [|J with 
G being the angle between the local direction and the nematic axis. It thus 
follows that in this case Tr an- = p/(l + F) and as announced, thus differs 
from the 3D equivalent Tr = p. The difference is due to the projection 
of the continuity condition Eq. onto a fluctuating surface M. 



Putting this ansatz into the minimization conditions we are left with the 
following results 



1-S 



2\2 



1 



2(/3ep) {(3u a - f (F„ - F x 



2{/3e) (1 - S 2 ) 2 (l + FY 

P i l TS) | 2p q max ^ 



2(/3e) 2 (l + F)(l-S 2 ) (3e27r(l + F)' 

Clearly the solution of the above equations describes a second order surface 
orientational transition characterized by S = and S ^ 0. The critical point 
of the transition is given as a solution of 

2(3ep = — ^ r-. (37) 



The transition involves only the variables (3, e and p, and can thus be achieved 
by either a variation in temperature, polymer stiffness or polymer surface 
coverage. 



3.6 Effective Elastic Constants 



The corresponding renormalization of the effective membrane elastic con- 
stants can be obtained on this level of approximations from 

S(Q) — S (Q)+PeY^c a Ql + 2\J2^Ql + -Y,^QlQl = 

a a P afi 
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a a P a ,p 

(38) 

The renormalized values of the elastic constants can be read off as 

7 — * 7 + 2-^<7 Q + ec Q 

K c — > ^ + 1^. (39) 

One observes first of all that the elastic modulus of the "dressed" membrane 
changes because of the elasticity of the embedded polymer, while the surface 
tension changes because of two distinct contributions: one stemming from 
the fact that polymer is embedded into the surface and thus scales linearly 
with A, the other one stemming from the coupling between the curvature of 
the surface and curvature of the embedded chain. 



If I now write 
(3S(Q) = Lz ]l (Ql + Ql)+lz x S(Ql-Ql) + 

HPK e (Ql + Ql) 2 + (V,, (Ql + Ql) + a ± S (<% - Q\ 

(40) 

thus defining the quantities a\\ and cr± as well as Sy and I obtain 

E,| = / 5 7 + 2/3ec|| + 2^2A + /3n a ^^^a|| 

S± = 2(3e c ± + 2 f2A + Ax a (i ^ F) W (41) 



while 



P 

(Tii = a_L 



2(1 + F) 

(2p Qrnax \ 
WeMLTT)- C± ) 

C± = 2(/5e) 2 (l + F)(l-5 2 )' (42) 

It is clear from Eq. [E] that depending on the orientational order parameter 
in the nematic phase the elastic energy will be destabilized if En — Y,±S < 0. 
This destabilization of the surface energy is due to the existence of an easy 
axis, in the direction perpendicular to the nematic axis, characterized by 
a much smaller surface energy than in the direction parallel to the nematic 
axis. This effect follows directly from the anysotropic packing of the polymer 
on the embedding surface. 
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3.7 The Tubulization Transition 



Renormalized elastic constants are now taken as a point of departure for 
the general analysis of the effective elastic energy. The recent quite general 
formalism developed by Radzihovsky and Toner || is taken as a lead. 



Instead of choosing the Monge parametrization for the surface r(x, y) = 
(x, y, ((x, y)) one can write down the elastic energy Eq. directly in terms 
of r(xi) with xi = x, x<i = y. Also in order to insure stability, in view of the 
discussion preceeding this section, appropriate fourth order terms have 
to be added to the expansion Eq. [H]. In the spirit of the Landau theory 
we assume that the fourth order displacement terms do not change sign as a 
function of the nematic order parameter, and thus assume them as constant. 



The total elastic energy, meaning the sum of Eq. [H] and the fourth order 
terms, can now be written in terms of the gradients of local displacements 
as i, 

F(rfa)) = \ J dx ± dx n K±± (^r^)) 2 + K nn (d 2 r(x t )) 2 + K n± (firfa)) (^r(^)) + 
+ t ± (d ± r(xi)) 2 + t n (d n r(xi)) 2 + 

+ ^u ±± (d±r(x i ).d_ L r(x i )) 2 + \u nn (d n r(xi).d n r(xi)) 2 + u n± (d n r(x i ).d ± r(x i )) 2 + 
+ ^v±± (d±r(xi).d ± r(xi)) 2 + v nl _ (d n r(xi)) 2 (d±r(xi)) 2 ] , (43) 

where we have chosen the nematic direction to be n. In order for the system 
to be stable the fourth order constants have to be positive u, v > 0. The 
second order constants can be obtained by comparison with Eq. [|(]. The 
effective elastic moduli are obtained as 

K ±± =K c + ^(a 2 n - a n a ± S + a\S 2 
P v 

K nn = K c + — [a 2 n + a n a ± S + a\S 2 



P 

K nL =K c + fj(oZ- a 2 ± S 2 ) , (44) 

and as can be easily deduced from the definitions Eq. f|2] are all positive. 
The effective surface tension can be cast into the form 

t n = S n + H±S 

t ± = E n - E X S, (45) 

where t± can be either positive or negative. The zero of t± does not coincide 
with the isotropic - nematic transition of the polymers, i.e. with S = 0. 
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This free energy is now treated in analogy with the usual 4 theories of 
the critical phenomena, identifying the tangent vectors t a = d a r(xi) as order 
parameters |J . In what follows we will delimit the analysis to the mean-field 
approximation where we take the ansatz tmf = (Ci x )Cn2/)> where the sym- 
metry breaking nematic axis has been oriented in the y (parallel) direction. 
The pref actors C_L, ( n are order parameters that measure the shrinkage of the 
membrane due to undulations fl(| . With this ansatz one obtains for the free 



energy 

T = \L ± .L n [ Un + t±tl + \ + v X x) Cl + lu nn (n + v n ±(n(l } ■ (46) 

Following Radzihovsky and Toner || by minimizing the above free energy 
one obtains the following phase diagram topology for (u±± + vj_±) u nn > v%±- 
There exists a crumpled phase with ( n = (j_ = for t n , t± > 0. As t± changes 

sign a tubular phase sets in with (± = and ( n = 0. This transition 

is second order. A transition to the flat phase, which can be expected on 
general grounds 0, never sets in as t n can not change sign. 



That £j_ is in deed bound to change sign somewhere within the polymer 

surface nematic phase is evident from the fact that Urns' >o^_l = S n ~ 

const. > while in the opposite limit linis ,it± = E n — T,± ~ —(1 — 

S 2 )~ < 0. The exact position of the change of sign does not coincide with 
the isotropic - nematic transition of the embedded polymer chain, i.e. with 
S = 0. Thus a sufficient ordering of the polymer chain has to be present 
before the tubulization transition of the whole embedding membrane can 
take place. 



A second order orientational transition of the embedded polymer chain thus 
drives an associated shape transition of the membrane, corresponding to a 
tubulization of the membrane in the direction perpendicular to the nematic 
axis, or parallel to the easy axis of the surface energy. The two transitions 
(i.e. the nematic to isotropic transition of the polymer and the symmetry 
breaking shape transition of the dressed membrane) do not, however, happen 
at the same value of the temperature (surface coverage or polymer stiffness). 



4 Discussion 



Strong adsorption of polymers onto soft supporting surfaces has different 
characteristics from ordinary polymer adsorption where the adsorption phe- 
nomenon results from a competition between adsorption energy and entrop- 
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ically driven tendency for less tightly bound configurations. Recent exper- 
imental work |3| shows that in the strong adsorption case the polymer 
chains are basically confined to lie on the adsorbing surface, without any 
dangling end or train configurations. In this case, with polymer chains ef- 
fectively embedded (onto) into the supporting surface, the main source of 
entropy is limited to either polymer conformational degrees of freedom along 
the embedding surface or to thermally driven fluctuations of the supporting 
surface. 



Because of this tight coupling of the polymer with the supporting surface, 
polymer statistics is strongly influenced by local surface configurations. Fur- 
thermore, if the polymer is non-ideal, i.e. self-interacting, the supporting 
surface fluctuations renormalize the interactions between polymer segments 
as has been amply demonstrated in a slightly different context by Goulian et 



al. fTT |. This effect depends on the nature of the polymer self-interactions 
as well as on the magnitude of the supporting surface fluctuations. We have 
investigated the properties and phenomena connected with an orientational 
ordering transition of surface-embedded polymers exhibiting an orientational 
(nematic) interaction potential and have shown explicitely the nature of the 
renormalization of the polymer self-interaction due to the supporting surface 
fluctuations. The conclusion reached in this connection is that for a suffi- 
ciently stiff membrane this effect tends to be small since it depends on the 
equilibrium fluctuations of the bare surface curvature. 



It is however not only true that thermally driven undulations of the sup- 
porting surface affect the polymer ordering transition, but this transition 
in its turn, manifestly modifies the properties of the compound (bare mem- 
brane + embedded polymers) membrane. The reverse effect, of polymer 
modified membrane properties, appears to be much more important both in 
qualitative and quantitative terms. Polymer orientational ordering spawns 
an associated symmetry breaking transition in the mean shape of the com- 
pound membrane that preferrs membrane configurations with most of the 
surface area in the direction of the easy (low energy) axis of the surface en- 
ergy. This is achieved through a tubulization transition where the long axis 
of the tubule is perpendicular to the nematic axis of the polymer ordering. 



Though it is understandable that nematic ordering of a surface embedded 
polymer should have some consequences on the properties of the "dressed" 
surface, the drastic change in the mean shape of the compound membrane 
does come as a surprise. 



One can foresee several effects that are beyond the presently formulated 
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mean-field solution to the tight adsorption model. First of all the coupling 
between local curvatures and the nematic order parameter of the polymer S is 
missing. It is intuitively plausible and in deed probable, that local curvature 
should affect local ordering of the polymer segments tending to concentrate 
alligned polymers in regions with large local curvature. This effect could 
also introduce additional ordered phases ]12| into the phase diagram of the 
"dressed" membrane. As is clear from comparison with the general analysis 
g] the ph ase diagram of the polymer embedding surface covers only one part 
of the phase diagram expected on general grounds. 



The tubulization transition, if indeed present in the systems with strong 
adsorption of polymers to soft flexible surfaces, would be of particular impor- 
tance for understanding the DNA - cationic lipid aggregation 0, [|, where 
the colloidal state of the aggregate is crucial for transfection of the DNA - 
lipid complex across the cellular membrane. 
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